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The two-dimensional Hubbard model exhibits superconductivity with d-wave symmetry even at 
half-filling in the presence of next-nearest neighbor hopping. Using plaquette cluster dynamical 
mean-field theory with a continuous-time quantum Monte Carlo impurity solver, we reveal the 
non-Fermi liquid character of the metallic phase in proximity to the superconducting state. Specif- 
ically, the low- frequency scattering rate for momenta near (7T, 0) varies non-monotonously at low 
temperatures, and the dc conductivity is T-linear at elevated temperatures with an upturn upon 
cooling. Evidence is provided that pairing fluctuations dominate the normal-conducting state even 
considerably above the superconducting transition temperature. 
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The main postulate of Landau's Fermi-liquid (FL) the- 
ory [l[ is the existence of well-defined fermionic quasi- 
particles with a lifetime which diverges more strongly 
than (£fc — Ep)^ 1 when their energy Ek approaches the 
Fermi energy Ep. The search for instabilities of the FL 
towards alternative metallic states which are not adia- 
batically connected to the non-interacting Fermi gas has 
developed into one of the central topics in correlated elec- 
tron physics. Perturbation theory of electron-electron in- 
teraction in dimensions higher than 1 reproduces the FL, 
implying that processes not contained therein are needed 
to destroy it [2j. We will present evidence in this letter 
that in the two-dimensional (2D) Hubbard model (HM), 
a generic model for interacting electrons on a lattice, pair- 
ing fluctuations are such processes which eventually lead 
to a breakdown of FL theory even for moderate inter- 
action strength in the proximity to a superconducting 
groundstate. 

The 2D HM has been studied extensively in the 
intermediate-to-strong coupling regime believed to be 
relevant for the physics of doped Mott insulators and 
cupratc superconductors. In the cuprates, deviations 
from FL behavior in the metallic regime were observed 
both in quantum oscillation and angle-resolved photoc- 
mission experiments 0, HI . Similarly, in the 2D HM the 
FL may become unstable towards antiferromagnetism 
(AF) [1, El or d-wave superconductivity (dSC) @-[loj]. 
Non-Fermi liquid (NFL) signatures and anisotropic pseu- 
dogaps have been identified for electron densities below 
half- filling Ill4l4|. As we discuss in this letter, a state 



with properties reminiscent of the strongly correlated 
state below half-filling already appears in the seemingly 
simpler weakly coupled half- filled system with finite next- 
nearest neighbor hopping t' . 

The possibility of a dSC instability at moderate cou- 
pling even at half-filling was previously discussed on the 



basis of renormalized mean-field 1 151 as well as cluster 
dynamical mean- field calculations [16| . A weak-coupling 
rcnormalization group analysis confirmed the existence 
of a dSC instability with a critical temperature which 
vanishes exponentially towards zero interaction strength 
. Functional renormalization group studies indicated 
the possibility of NFL normal state properties, in prox- 
imity to van Hove singularities 
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Here we study the half-filled 2D HM in the paramag- 
netic phase at finite but low temperatures. Using pla- 
quette cluster dynamical mean- field theory (CDMFT), 
we confirm the existence of dSC at moderate coupling. 
We investigate the normal-state properties above the 
dSC critical temperature and reveal an unusual metallic 
state characterized by an increasing instead of decreas- 
ing single-particle scattering rate upon lowering the tem- 
perature T. Moreover, the dc resistivity is T-linear at 
elevated temperatures and shows an upturn at lower T. 
This evidence either implies that dSC emerges from an 
unusual metallic state with NFL character or that super- 
conducting fluctuations themselves severely modify the 
metallic state above T c . 



The 2D U-t-t' Hubbard model reads 



H 



(1) 



where ct (c fe ) creates (annihilates) an electron in 
a Bloch state with momentum k and dispersion = 
—2t (cos k x + cos k y ) — it' cos k x cos k y , rii^ G is the density 
operator for site i and spin a = f, |, U > is the local 
Coulomb repulsion strength, and \x the chemical poten- 
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FIG. 1. Overview of paramagnetic phases obtained for the 
half-filled 2D U-t-t' Hubbard model within 2x2 CDMFT at 
T/t — 1/100. Data points show metallic (blue squares), dSC 
(red circles) and Mott insulator (black diamonds) states. For 
the characterization of phases see Fig. [2] The shaded areas 
are a guide to the eye. 



FIG. 2. The superconducting order parameter ^(tt,o) (left 
axis) and the double occupancy D (right axis) as a function 
of U/t. For U/t > 6.5 only the Mott-insulating solution to the 
CDMFT equations is shown, but note that a metallic solution 
also exists in a hysteresis region of U values. 



tial. The CDMFT self-consistency equations @, are 
GQu n ) = J2 + A 1 _ ~ *(*)) 1 ' ( 2 ) 

k 

^o 1 ^) = G _1 (icj„) - 27(iw n ). (3) 

For the N c = 2 x 2 plaquette, the hopping ma- 
trix t(k) is defined via its matrix elements tij (k) = 



(k+k)-(Xi-Xj 



e fc+fc' wnere -X* an d Xj are the 



position vectors of cluster sites i and j, fe is in the re- 
duced Brillouin zone, and the cluster momentum takes 
the values k = (0,0), (n, 0), (0,7r), and (tt, 7r). All quan- 
tities, i.e. i, the coarse-grained cluster Green function 
G, the corresponding Weiss field Go and the cluster self- 
energy £ are N c x 7V C matrices, and 1 is the unit matrix. 
In the following we consider only paramagnetic solutions 
and the spin index is therefore suppressed. The self- 
consistency cycle is closed by solving the impurity prob- 
lem, i.e. by calculating, for a given self-energy and the 
corresponding Weiss field, a new cluster Green function 
matrix Gij(r) = — (7^Cj(r)ct (0))s cff , where S c g denotes 
the effective action of the auxiliary Anderson impurity 
model. 

The CDMFT formalism is generalized to the su- 
perconducting state by introducing the spinors \Pk = 
(Cfe-j-, c_ fe i ). The corresponding Nambu Green function 
matrix in momentum space reads 



G fc (r) = (%(tM) 



G fe ,t(r) F fc (r) 



(4) 

Within plaquette CDMFT the cluster Green function 
matrix is diagonal in cluster momentum space, and the 



relevant anomalous Green functions for ri-wave order 
parameter symmetry are F^^(t) = —F(q^(t), while 
F( ,o) an d F(n,ir) both vanish. If a self-consistent dSC so- 
lution is obtained, the order parameter is = |-Ffc(r = 
+ )l = l(c fct (0+) C _ fci )|forfc = (7r,0). 

The impurity problem is solved by numerically exact 
continuous-time quantum Monte Carlo (QMC) simula- 
tions based on the expansion of the effective action in the 
impurity-bath hybridizationfiij generalized to supercon- 
ducting states [25[. The hybridization expansion QMC 
allows to reach temperatures sufficiently low to detect 
dSC instabilities within plaquette CDMFT and to trace 
pairing fluctuations to temperatures well above T c . 

We start with an overview of the plaquette CDMFT 
phases obtained for T/t = 0.01 in Fig. [TJ Four distinct 
low-temperature phases are indicated: a metal at very 
weak coupling, a superconducting state with ri-wave or- 
der parameter symmetry at moderate coupling, a metal 
with NFL properties, and a Mott insulator at strong cou- 
pling. Below we focus specifically on \t' /t\ = 0.3 and an- 
alyze the distinctive character of the individual phases. 

The average double occupancy D = N^ 1 J2iti ( n it n il) 
decreases with increasing U, until the discontinuous Mott 



transition [2 6H2 8j is reached where a sharp drop of D at 
U c between 6.25 and 6.5 occurs (see Fig. [5]). A non-zero 
superconducting order parameter |Z\( w )| is obtained in 
a converged dSC solution of the CDMFT equations in 
the range 1 < U/t < 5. In the metallic phases the dSC 
order parameter vanishes. For even lower temperatures 
T < t/100 the dSC region in the phase diagram is ex- 
pected to grow and most likely to extend to arbitrarily 
weak interaction strengths in the zero-temperature limit 
- in conjunction with asymptotically exact weak-coupling 
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FIG. 3. Imaginary part of the self-energy at k = (0, 0) and 
(it, 0) at the Matsubara frequencies iu„ = i(2n + 1)ttT for U/t 
= 4, \t' /t\ — 0.3, and different temperatures. Even though 
a dSC state is found for these parameters at T/t = 1/100 
we show here the self-energies for a solution of the CDMFT 
equations where only the diagonal Green functions are non- 
zero in Eq. @. 



rcnormalizaton group results [17 1 . 

In Fig.[3]we show the imaginary part of the self-energy 
Im£k(iw n ) for the cluster momenta k = (0, 0) and (tt, 0) 
for U = At in the metallic phase in proximity to super- 
conductivity. At k = (0, 0), i.e. far away from the Fermi 
surface, the absolute value of ImS grows upon lowering 
the temperature below T/t = 1/50; dSC sets in below T/t 
ss 1/90. The same behavior is also observed for stronger 
coupling U/t = 5 and for t' = at U/t = 4. 

The self-energy at (tt, 0), the closest point to the Fermi 
surface resolved in plaquette CDMFT, instead is remark- 
ably different. In a 2D Fermi liquid the generic low- 
frequency form of the self-energy [29J ImU '(uj) oc uj 2 lnw + 
T 2 lnT implies that also ImI7(iu;„) has a negative slope 
for uj n — > 0, with a zero- frequency offset which van- 
ishes for T — >• 0. In contrast, the self-energy at fc = 
(tt, 0), shown in the lower panel of Fig. [31 develops a 
low-frequency hump feature with a change in slope upon 
lowering T, indicative of an increased scattering rate near 
the Fermi energy. This hump feature is distinct from 
Fermi liquid behavior; it emerges uniquely in the param- 
eter regime in which dSC occurs. For example, if t' = 
or U/t > 5 the dSC instability and the NFL hump feature 
in the self-energy at (tt, 0) are both absent. Within our 
numerical approach it is, however, difficult to trace this 
low-energy NFL anomaly to even lower T or smaller U 
values. At k ~ (tt,tt), Im E(iui n ) is almost independent 
of temperature and vanishes linearly at low frequencies. 

These features for the normal-state self-energy suggest 
unusual transport properties, too. We therefore calcu- 
late the dc resistivity by performing a convolution on 



the imaginary frequency axis without vertex corrections 
[30j . The longitudinal dc conductivity is obtained ap- 
proximately via 



Ode 



T 2 TTh 



A xx (t = 1/2T) 



(5) 



from the long- wavelength limit of the current-current cor- 
relation function 

a xx (t) = ^j2^ e ^ 2e ~ iUmTG ^+^ G ^ iu M 

nmh 

for bosonic Matsubara frequencies u> m = 2rmrT. N 
is the number of fc points in the first Brillouin zone, 
G h (iL> n ) = (iuj n + fi - e h - Zfeflo^)) -1 , and S k = 
\ Ylij e lk '( Xi ~ x i^ Sij is an interpolated lattice self-en ergy 
obtained from the cluster impurity self-energy Eij [221 ]. 
The basis for the approximate Eq. (JS[) is the relation be- 
tween the current-current correlation function at imagi- 
nary time r and the imaginary part of the current-current 
correlation function at real frequency u>, 



A xx (t) = 



did cue UJT ImA xx (uj) 



Loo tt 1- e~»/ T 



(7) 



Within CDMFT A xx (t) can be directly obtained from 
Eq. ([6]), but to determine Ode it is necessary to com- 
pute lim w _ ) .o lmA xx (ui)/u! [3l[. However, T(\-e-"/ T ) l& 
an even function, strongly peaked at lo/T = and decays 
exponentially in u/T. Therefore, T^A^t = 1/2T) w 
lim w _j.o ImA xx (u>) /ui asymptotically for T — > 0; hence Eq. 
(0 follows. 

The resulting T-dependent dc resistivity pdc = &dc f° r 
selected values of U ft is shown in Fig. 0] Remarkably, 
the resistivity is T-linear at elevated temperatures. A 
deviation from the T-linear resistivity with an upturn at 
lower temperatures sets in at T/t « 0.02 for all U val- 
ues. Therefore, adding to the anomalous behavior of the 
antinodal (tt, 0) self-energy, also the T-linear resistivity 
and the low-temperature upturn are different from Fermi 
liquid behavior. Pdc(T) develops its upturn at approxi- 
mately the same temperature at which the hump feature 
in the antinodal self-energy appears. 

This hump bears resemblance with a similar feature 
observed in the self-energy on the real-frequency axis in 
the fRG results of Ref. [l9|. The spike structure in Ref. [3 
is especially pronounced near points on the Fermi surface 
which are connected by the umklapp scattering wave vec- 
tor Q = (tt, 7t). The origin of the spike feature in the fRG 
study is therefore attributed to enhanced umklapp scat- 
tering due to the emerging antiferromagnetic spin fluc- 
tuations. Similarly, also the NFL hump feature in our 
study suggests an explanation in terms of an emergent 
scattering mechanism. Since the (tt, 0) self-energy devel- 
ops its NFL structure uniquely for parameter values close 
to the dSC state, it is suggestive that they are caused by 
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FIG. 4. T-dependence of the dc resistivity pdc (without ver- 
tex corrections) for selected values of U/t. At high T the 
resistivity is nearly linear in T, as indicated by the dashed 
lines (piin) which connect the origin and the data points at 
the respective highest temperatures. Deviations from this T- 
linear behavior set in below T/t fs 0.02. Note that each curve 
stops at the smallest temperature at which the solution is not 
superconducting. Inset: Difference between pdc and pn n . 



superconducting fluctuations, similarly as in the attrac- 
tive Hubbard model in proximity to the superconducting 
state (3^ |. Additional insight regarding the importance 
of superconducting fluctuations above T c may be gained 
from the probabilities of relevant plaquette cigenstates 
in the thermal ensemble 



33 1 . We have verified that a 



spin singlet two-hole state with zero total momentum, 
i.e. a "Cooper pair state" with respect to the half-filled 
plaquette, is represented with enhanced probability not 
only in the dSC state at low temperature, but already 
significantly above T c . 

Despite the limitations imposed by neglecting current- 
vertex corrections and by the restricted momentum-space 
resolution, the T-linear behavior of the dc resistivity at 
elevated temperatures and its upturn at low tempera- 
tures mark nevertheless clear deviations from Fermi liq- 
uid physics and bear resemblance to resistivity data in 
cuprates [34j and in other families of high-temperature 
superconductors as well j35|. Specifically superconduct- 
ing phase fluctuations were recently verified in thermo- 
dynamic measurements well above T c and discussed in 
connection with the unusual properties of the metallic 
phase in the cuprates [36| . The resistivity upturn we find 
for the half-filled Hubbard model may therefore be tied 
to the superconducting fluctuations in the vicinity of the 
dSC phase. We note that the resistivity upturn occurs 
also for U/t = 5, but at a somewhat lower temperature 
than for the smaller U/t values. This is in accordance 
with Fig. [5J which shows that if superconductivity oc- 
curs for U/t = 5, the critical temperature is smaller than 



T/t = 0.01. 

Since we focused here on finite-temperature results 
for the 2D model, long-range AF order was naturally 
not considered, while dSC can emerge via a Kostcrlitz- 
Thoulcss type transition |37j. At T = 0, AF is strongest 
for t' = with the lower critical coupling strength U^ F 
= [1,0]. For \t'\ > 0, long-range AF order is only ob- 
tained beyond a critical coupling U^ F (t') > p|. The 
RG results in Ref. [I?] imply the existence of a finite dSC 
range for \t'\ > and weak to moderate U at T = 0. At 
U^ F we expect a transition to a ground state with AF 
order, possibly in coexistence with dSC [la. Ilo. l38|. 

Our findings show that even for moderate correlations 
the physics of the half-filled 2D Hubbard model with t' ^ 
already contains the superconducting and the unusual 
metallic NFL phases encountered in previous studies of 
the doped system at stronger coupling. It seems likely 
that the strong-coupling dSC and NFL phases at finite 
hole doping are continuously connected to their counter- 
parts at weak coupling and half-filling for finite t' — a 
hypothesis to be elaborated on in future work. Finally 
wc note that our findings of NFL signatures also bear 
resemblance to the single-site DMFT results for the nor- 
mal state of the attractive Hubbard model close above 
the transition to s-wave superconductivity (39[. 
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